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Massive stars (with mass > 8M Q ) are fundamental to the evolution of galaxies, because 
they produce heavy elements, inject energy into the interstellar medium, and possibly regulate 
the star formation rate. The individual star formation time, t*f, determines the accretion 
rate of the star; the value of the former quantity is currently uncertain by many orders of 
magnitude 1 - 23 - 4 ' 5 ' 6 , leading to other astrophysical questions. For example, the variation oft,/ 

O 

with stellar mass dictates whether massive stars can form simultaneously with low-mass stars 

(N . 

O . in clusters. Here we show that U } is determined by conditions in the star's natal cloud, and 

is typically ~ I0 5 yr. The corresponding mass accretion rate depends on the pressure within 
the cloud — which we relate to the gas surface density — and on both the instantaneous and 
final stellar masses. Characteristic accretion rates are sufficient to overcome radiation pressure 
from — 100 M protostars, while simultaneously driving intense bipolar gas outflows. The weak 

X. 

dependence of <*/ on the final mass of the star allows high- and low-mass star formation to 
occur nearly simultaneously in clusters. 

Massive stars form in dense molecular clumps inside molecular clouds 7 . These regions are highly turbu- 
lent and are in approximate virial equilibrium, with comparable values of the gravitational energy and the 
kinetic energy. Observed star forming regions have a virial parameter a V i r = 5(a 2 )R c \/GM c \ of order unity 8 , 
where (a 2 ) is the mass-averaged one-dimensional velocity dispersion, R c i the radius, and M c \ the mass of the 
clump. The regions of high-mass star formation studied in ref. 7 are characterized by masses M c \ ~ 3800 M Q 
and radii R c \ ~ 0.5 pc. The corresponding mean column density and visual extinction are E ~ 1 g cm -2 and 
Ay — (Nn/2 x 10 21 cm -2 ) mag = 214 £ mag, with a dispersion of a factor of a few. These column densities 
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are far greater than those associated with regions of low-mass star formation, but are comparable to the 
central surface density in the Orion Nebula Cluster 9 , which is about 1 gcm~ 2 . These regions have very high 
mean pressures P: 



where k B is Boltzmann's constant and / gas is the fraction of the cloud's mass that is in gas, as opposed to 



Molecular clouds are observed to be inhomogeneous on a wide range of scales, and numerical simulations 
show that this is a natural outcome of supersonic turbulence 11 . Dense, self-gravitating inhomogeneities that 
are destined to become stars are termed cores. We assume that the distribution of core masses determines 
the resulting distribution of stellar masses (the initial mass function) and we take this distribution as given. 
Because of the disruptive effects of protostellar outflows, only a fraction e coro of the total core mass M core 
ends up in the star: m*/ = e C ore^corc, where to*/ is the final stellar mass (for binary and other multiple 
star systems, m*/ is the total mass of stars in the system). To estimate e C orc for massive stars, we assume 
that their outflows are scaled versions of the outflows from low-mass stars, which is qualitatively consistent 
with observation 12 . For low-mass stars, protostellar outflows typically eject 25-75% of the core mass 13 . A 
similar analysis applied to high-mass stars yields comparable results (J.C.T. & C.F.M., in prep.), so we 
adopt e coro = 0.5 as a typical value. 

Our fundamental assumption is that star-forming clumps and the cores embedded within them arc 
each part of a self-similar, self-gravitating turbulent structure that is virialized — that is, in approximate 
hydrostatic equilibrium — on all scales above the thermal Jeans mass. The density and pressure are then 
power laws in radius (for the pressure, P oc r~ kp ), so that the cloud is a polytrope with P oc p lp . In 
hydrostatic equilibrium 14 , k p = 2/(2 - j p ) and k P = lp k p = 2 7p /(2 - 7p ). Let c = {P/p) 1/2 oc A 1 -^^ 2 -^ 
be the effective sound speed; if the pressure is dominated by turbulent motions, then c is proportional to the 
line width. Molecular clouds and cloud cores satisfy a line width-size relation in which c increases outwards 15 , 
corresponding to 7 P < 1. The equation of hydrostatic equilibrium gives M — kpc 2 r/G for the mass inside 
r, and p = Ac 2 j (2^Gr 2 ) for the density at r, where A — 7 P (4 — 37 P )/(2 — 7 P ) 2 . ft is then immediately 
possible to determine the properties of a core in terms of the pressure at its surface and the mass of the 
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stars. They have power-law density profiles 



p cur 



k », with k p ~ 1.5 ±0.5. 



star that will form within it (see Fig. 1). The radius of a core is then r s = 0.074(to*//30 Mq) 1 /^ -1 / 2 pc; 
recall that the typical clump observed in ref. 7 has a radius of 0.5 pc, which is set by the angular resolution 
of those observations. The mean density in a core is nn = 1.0 x lO^m^/SOMo)- 1 / 2 ^/ 2 cm , and the 
r.m.s. velocity dispersion is 1.65(to*//30 M Q ) 1 / 4 S 1 / 4 krns" 1 . The thermal sound speed at temperature T is 
0.3(T/30Ky/ 2 kms-\ and so a high-mass core is dominated by supersonic turbulence and should be very 
clumpy. 

We now consider the timescale for a star to form in such a core. On dimensional grounds, we expect 
the protostellar accretion rate to be 16 : 

m* = 0* — , (2) 



so long as radiation pressure is not important. Here to* is the instantaneous stellar mass, tg = [3n/ (32G/})] 1//2 
is the free-fall time and 0* is a dimensionless constant of order unity. This equation could be violated in 
the sense that 0*^1 only in the unlikely case that the star forms from a coherent spherical implosion; 
if the star formation is triggered by an external increase in pressure, </>* could be increased somewhat, but 
deviations from spherical symmetry in the triggering impulse and in the protostellar core will generally 
prevent <j)„ from becoming too large. It could be violated in the opposite sense that 0* <C 1 if the core is 
magnetically dominated, so that collapse could not begin until the magnetic field diffused out of the core. 
However, magnetic fields are not observed to be dominant in molecular cores 17 . Equation (2) implies that 
the accretion rate, and thus the star formation time t*/ oc to*/to*, depends weakly on the properties of the 
ambient medium, to* oc p 1 / 2 . 

We now show that if the collapse is spherical and self-similar, then </>* is quite close to unity provided 
that radiation pressure does not disrupt the flow. Although we assume the collapse is spherical far from the 
star, it will naturally proceed via a disk close to the star owing to the angular momentum of the accreting 
material; we assume this does not limit the flow of matter onto the star as otherwise the disk would become 
very massive and gravitationally unstable. Shielding by the disk reduces the importance of radiation pressure 
on the accretion flow, and allows the formation of massive stars provided the accretion rate is sufficiently 
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large 18 . Under the assumption of a polytropic structure in hydrostatic equilibrium, equation (2) implies: 

4 c 3 

TO* = 0*£ C orc 7^ fc P-4 1/2 -— (3) 

7TVO 

where the value of p entering tg is evaluated at the radius in the initial cloud that just encloses the gas that 
goes into the star when its mass is to*. For the isothermal case (j p — 1) and for e cor o = 1, this reduces to the 
classic result of Shu 19 with 0* = 0.9757rV3/8 = 0.663. McLaughlin & Pudritz 2 show that the accretion rate 
to* is proportional to i 3 ~ 37p , so that for -f p < 1, as observed, the accretion rate accelerates 2 ' 4 . As discussed 
in ref. 2, termination of the accretion breaks the self-similarity once the expansion wave encloses to*/, which 
occurs at t ~ 0.45i*/. Thereafter, the relation to* oc i 3 ~ 3 Tp becomes approximate, but the approximation 
should be reasonable good 2 . The star- formation time is then given by: 

U f - ^^te- (4) 

Evaluating (/>* (see Fig. 1), we conclude that to* ~ m*/tg to within a factor 1.5 for spherical cores in which 
the effective sound speed increases outward. 

We can express the accretion rate in terms of the mean pressure of the clump in which the star forms, 
P (eq. 1), and the final mass of the star, to*/, 

to* = 4.75 x 10- 4 eV r 4 e(/ gas ^avi r ) 3/8 (^J^) 3 ^ 374 (^J Moyr" 1 , (5) 

where <pp = P s /P is the ratio of the core's surface pressure to the mean pressure in the clump and j = 
3(2 — k p )/[2(3 — k p )]; for k p = 3/2, j = 1/2. This typical accretion rate is large because the core is embedded 
in a high-pressure environment, so that in hydrostatic equilibrium it has a high density and short free- 
fall time. Because massive cores are turbulent and clumpy, we expect the accretion rate to exhibit large 
fluctuations. Whereas dumpiness is an attribute of massive star-forming regions in our model, it is not a 
prerequisite, as suggested in ref. 3. The star formation time is: 

1/4 



Uf = i.26 x wxo\iHf s ^ P ^r 3/8 (^) £- 3 /v. 



(6) 



The weak dependence of the star formation time on the final stellar mass means that low-mass and high- 
mass stars can form approximately at the same time. Furthermore, i*/ is small compared to the estimated 
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timescale ~ 10 6 yr for cluster formation 20 . With the fiducial values of the parameters, a 100 M Q star forming 
in a clump with T, = 1 gem -2 has a final accretion rate of 1.1 x 10~ 3 M Q yr -1 and a star formation time of 
1.8 x 10 5 yr. 

A necessary condition for massive star formation is that the ram pressure associated with accretion 
exceed the radiation pressure at the point where the dust in the infalling gas is destroyed 21 . For a 100 M 
star, this requires 18 m»^6x 10~ 4 M Q yr -1 , which is satisfied for S ^ 0.5 g cm~ 2 . Once the core has formed 
a star, the star can continue to grow by Bondi-Hoyle accretion 22 provided m* <; 10 M Q so that radiation 
pressure does not prevent the focussing of gas streamlines in the wake of the star. The Bondi-Hoyle accretion 
rate m* = 2.2 x 10- 7 (M cl /4 x 10 3 M Q )- 5 / 4 E 3 / 4 (m*/10 M Q ) 2 M Q yr" 1 (CFM & JCT, in prep.) is so low 
that this process does not significantly alter the stellar mass under the conditions observed in the clumps of 
ref. 7. 

Direct comparison of our results with observation is difficult because the actual masses of massive 
protostars are poorly determined. Our approach is to predict the properties of some well-studied massive 
protostars in terms of their bolometric luminosities. The bolometric luminosity Lboi has contributions from 
main-sequence nuclear burning L ms , deuterium burning Lp, and accretion L acc . The accretion luminosity 
L acc = /accGm»m»/r„ where / acc is a factor of order unity accounting for energy radiated by an accretion 
disk, advected into the star or converted into kinetic energy of outflows, and where the stellar radius r* may 
depend sensitively on the accretion rate m*. Massive stars join the main sequence during their accretion 
phase at a mass that also depends on the accretion rate 23 . To treat accelerating accretion rates, we have 
developed a simple model for protostellar evolution based on that of Nakano et al. 6,24 The model accounts 
for the total energy of the protostar as it accretes and dissociates matter and, if the central temperature 
T c <; 10 6 K, burns deuterium. We have modified this model to include additional processes, such as deuterium 
shell burning, and we have calibrated these modifications against the more detailed calculations of Palla & 
Stahler 25 ' 23 . 

Our model allows us to make predictions for the masses and accretion rates of embedded protostars 
that are thought to power hot molecular cores (CFM & JCT, in prep.). Figure 2 compares our theoretical 
tracks with the observed bolometric luminosities of several sources. We find that uncertainties in the value 
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of the pressure create only small uncertainties in to* for iboi k, few x 10 4 L Q . 

The infrared and submillimeter spectra of accreting protostars and their surrounding envelopes have 
been modelled in ref. 5., modelling the same sources shown in Fig. 2. We note that uncertainties in the 
structure of the gas envelope and the possible contributions from additional surrounding gas cores or diffuse 
gas will affect the observed spectrum. Comparing results, our inferred stellar masses are similar, but our 
accretion rates are systematically smaller by factors of ~ 2 — 5. The modelled 5 high accretion rates of 
~ 10 _3 M Q yr _1 for stars with to* <~ 10 M Q would be difficult to achieve unless the pressure was increased 
substantially; for example, if the stars are destined to reach to*/ ~ 30 M , pressure increases of a factor 
~ 40 are required. 
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Fig. 1. Variation of model parameters and results with k p . Let P s = <fipP be the surface 
pressure of a core. The properties at the surface of the core are given by r s = [AGm 2 ^ / (2irk P el OIC P s )] 1 ^ 4 , 
p s = [Ak 2 P e 2 COIC P 3 1 (27rG 3 ml f )} 1/4 , and c s = [2i:G 3 ml f P s /(Ak 2 P e 2 OIC )} 1/8 . We anticipate that the overall star 
formation efficiency in the clump will be relatively high, so in equation (1) we adopt / gas = 2/3 as a fiducial 
value. We estimate cf)p = P s /P ~ 2 (C.F.M. & J.C.T., in prep.), and we set a v ir equal to unity 7 . We take k p = 
1.5 as a typical value, corresponding to 7p = 2/3, kp = 1 and A = 3/4. Following equation (3), we evaluate 
0* using the results of McLaughlin & Pudritz 2 , </>* = tt\/3[(2 - 7p ) 2 (4 - 3 7p )( 7 - 6 ^/ 2 8( 3 ^- 5 )/ 2 m ] 1/(4 - 37p) , 
where too is tabulated in their paper. Over the entire range of 7p and k p relevant to molecular clouds 
(0 < 7 P < 1, 1 < k p < 2), (f>* ~ 1.62 — 0.48fcp to within about 1%. The star formation time decreases from 
3.5iff to 1.5iff as 7p varies from to 1. The variation of t*f and m*f relative to the k p = 1.5 case is also 
shown. Note that the singular polytropic model in hydrostatic equilibrium breaks down for k p = 1, 7p = 0, 
since then the pressure gradient vanishes (kp = 0). 

Fig. 2. Derived properties of nearby massive protostars. Solid lines show the predicted evolution 
in luminosity of protostars (including their accretion disks, but allowing for the powering of protostcllar 
outflows so that / acc = 0.5; J.C.T. & C.F.M. , in prep.) of final mass 7.5, 15, 30, 60 and 120 M Q accreting 
from cores with k p = 1.5 embedded in a S — 1 g cm -2 clump, typical of Galactic regions 7 . The luminosity step 
occurring at around 5 to 7 M Q , depending on the model, corresponds to the onset of deuterium shell burning, 
which swells the protostellar radius by a factor of about two and thus reduces the accretion luminosity by 
the same factor. The dashed and long dashed lines show a 30 M© star forming in a clump with mean 
pressure ten times smaller and larger than the fiducial value, respectively. The dotted line shows the zero age 
main sequence luminosity from ref. 26. Four observed hot molecular cores are shown: G34.24+0.13MM 27 , 
IRAS 23 3 85+60 53 28 , Orion Hot Core 29 and W3(H 2 0)— the Turner-Welch object 30 . The vertical error bar 
illustrates the uncertainty in the bolometric luminosity. The horizontal error bar shows the corresponding 
range of allowed values of to* for the S = lgcm~ 2 models. These values and the constraints on to* are 
listed here for each source. 
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